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ABSTRACT 

We discuss alternative interpretations of the differences in the Spectral Energy Distri- 
butions (SEDs) of BL Lacs found in complete Radio or X-ray surveys. 
A large body of observations in different bands suggests that the SEDs of BL Lac 
objects appearing in X-ray surveys differ from those appearing in radio surveys mainly 
in having a (synchrotron) spectral cut-off (or break) at much higher frequency. 
In order to explain the different properties of radio and X-ray selected BL Lacs Giommi 
and Padovani proposed a model based on a common radio luminosity function. At 
each radio luminosity, objects with high frequency spectral cut-offs are assumed to 
be a minority. Nevertheless they dominate the X-ray selected population due to the 
larger X-ray-to-radio-flux ratio. An alternative model explored here (reminiscent of 
the orientation models previously proposed) is that the X-ray luminosity function is 
"primary" and that at each X-ray luminosity a minority of objects has larger radio-to- 
X-ray flux ratio. 

The predictions of the two scenarios, computed via a Montecarlo technique, are com- 
pared with the observed properties of BL Lacs in the two samples extracted respectively 
from the 1 Jy radio survey and the Einstein Slew Survey. We show that both models 
can explain a number but not all the observed features. 

We then propose a completely new approach, based on the idea that the physical pa- 
rameter which governs the shape of the SEDs, is (or is associated with) the bolometric 
luminosity. Assuming an empirical relation between spectral shape and luminosity we 
show that the observational properties of the two surveys can be reproduced at least 
with the same accuracy as the two previous models. 

Key words: galaxies: jets, luminosity function - BL Lacertae objects: general - ra- 
diative mechanisms: non-thermal - surveys - methods: statistical 



1 INTRODUCTION 

Among Active Galactic Nuclei, Blazars show extreme lumi- 
nosity and variability. The observed properties are best in- 
terpreted in the frame of the relativistic jet model, proposed 
almost twenty years ago by Blandford & Rees (1978), as due 
to relativistic "beaming", i.e. the effects of the relativistic 
motion of the emitting plasma on the observed radiation. In 
particular BL Lac objects, which are considered part of this 
class, are characterized by their almost featureless continua. 

BL Lacs have been almost exclusively discovered 
through radio or X-ray surveys. However the properties of 
objects selected in the two spectral bands are systemati- 
cally different, posing a question as to whether there are two 



"types" of BL Lacs. The first difference to be recognized and 
perhaps still the most striking is the shape of the SED. The 
differences show up using broad band spectral indices and 
color-color diagrams, e.g. aao vs aox^ (Stocke et al. 1985; 
Maraschi etal. 1995; Sambruna, Maraschi & Urry 1996). 
Others include optical polarization, variability, presence of 
(weak) emission lines, radio luminosity, core-dominance, all 
of which are less conspicuous in X-ray selected objects (e.g. 
KoUgaard et al. 1992; Perlman & Stocke 1993; Jannuzi et al. 



* oi_2 = — log(^l/^2)/ log{i^i/i'2); radio fluxes arc taken at 5 
GHz, optical at 5500 A, X-ray at 1 keV 
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1994; KoUgaard et al. 1996; see KoUgaard 1994 for a recent 
review) . 

On the basis of the X-ray to radio flux ratio (or anx 
value) we can define an objective criterion (independent of 
the selection band) separating two (putative) classes of ob- 
jects: we define XBLs the objects with log(_Fikcv/^5GHz) > 
—5.5 and RBLs the objects with a ratio smaller than this 
dividing value, where the fiuxes are monochromatic and ex- 
pressed in the same units (see also Wurtz 1994, Giommi & 
Padovani 1994, for an analogous definition). As can be seen 
from Fig. la XBLs are found mostly but not exclusively in 
X-ray surveys and the same is true for RBLs with respect 
to radio surveys. 

The spread in spectral shapes was originally attributed 
to orientation effects, associated with different widths of the 
beaming cones of the radio and X-ray radiation emitted by 
a relativistic jet (Stocke etal. 1985; Maraschi etal. 1986; 
Celotti et al. 1993) . The idea came from the observational 
evidence that BL Lacs discovered in radio and X-ray surveys 
actually show similar X-ray luminosities while the radio lu- 
minosities typically differ by two-three orders of magnitude 
(e.g. Fig. lb). This could be accounted for if X-ray radiation 
had a wider beaming cone than radio emission: observers 
would see similar X-ray luminosities over a wide range of 
angles while the accompanying radio luminosity would be 
high for a small fraction of objects seen at very small angles 
and strongly dimmed for the majority, observed at larger 
angles (e.g. Fig. 5 in Celotti etal. 1993). Consequently the 
number density ratio between the two "fiavours" would be 
determined by the associated solid angles. Since the X-ray 
emission is largely isotropic, X-ray surveys are not biased 
against any of the two classes of objects, and can give the 
correct number ratio. The different beaming affecting the 
various bands could be due to an accelerating (Ghisellini 
& Maraschi 1989) or an increasingly coUimated jet (Celotti 
etal 1993). 

Substantial progress in the data on the SEDs of BL Lac 
objects has been obtained in recent years, in terms both 
of sensitivity and statistics, especially in the X-ray band 
(Giommi, Ansari & Micol 1995; Comastri, Molendi & Ghis- 
ellini 1995; Perlman etal. 1996a,b; Urry et al. 1996; Sam- 
bruna et al. 1996) . In terms of the power emitted per decade 
XBLs display a continuous rise up to the UV and in extreme 
cases the soft X-ray band, while RBLs are characterized by 
a spectral turn over in the IR domain. In XBLs the X-ray 
emission is dominated by a soft spectral component which 
extrapolates continuously to lower frequencies. In RBLs the 
X-ray emission is dominated by a separate harder compo- 
nent which at least in some cases extends to the 7-ray do- 
main (von Montigny et al. 1995; for X-ray spectra see the 
results by Padovani & Giommi (1996) and Lamer, Brunner 
& Staubert 1996). A hard GeV to TeV component is present 
also in XBLs though it does not show up in medium energy 
X-rays. Sambruna et al. (1996) showed that it is difficult to 
model the detailed transition from an XBL to an RBL like in 
terms of orientation only and suggested rather a continuous 
change in the physical parameters of the jet. 

The first model developed along this line to explain the 
"statistics" of XBLs and RBLs is the "radio luminosity -f 
different energy cutoff" scenario by Giommi and Padovani 
(Giommi & Padovani, 1994; Padovani & Giommi 1995; here- 
after we refer to it as the "radio leading" scenario by G&P) 



alternative in many ways to the commonly accepted "dif- 
ferent viewing angle" model sketched above. They propose 
that a single luminosity function in the radio band describes 
the full BL Lac population. For each radio luminosity X-ray 
bright BL Lacs (i.e. XBLs) are intrinsically a minority de- 
scribed by a fixed (luminosity independent) distribution of 
X-ray to radio flux ratios^ In X-ray surveys however selec- 
tion effects substantially enhance the XBL fraction. Accord- 
ing to this approach the intrinsic fraction of the two types 
of BL Lac (XBL vs. RBL) would be objectively reflected 
in radio surveys. With these hypotheses G&P were able to 
reproduce the observed X-ray counts, luminosity functions 
and the distribution of BL Lacs in the qro — Qox plane. 

Because of the implications of these issues for physical 
models of relativistic jets and the understanding of the phys- 
ical conditions within the emission region, we decided to ex- 
plore more thoroughly the fundamental hypothesis that BL 
Lacs are a single class, whose SEDs are characterized by dif- 
ferent physical parameters, and test any prediction against 
the observations now available. 

In addition to the radio leading model we consider a first 
"symmetric" alternative, namely that the X-ray luminosity 
function basically represents the whole BL Lac population. 
In this case X-ray surveys would give objective results re- 
garding the intrinsic abundance of XBLs and RBLs. This 
scenario, which we refer to as the "X-ray leading" model, 
can be considered as an evolution of the "different viewing 
angle" scenario, where the X-ray luminosity was the basic 
property and the population ratio reflected the ratio of solid 
angles. 

The predictions of both models, derived by a Monte- 
carlo technique, are compared in detail with the observa- 
tions now available. A main feature of the samples, that is 
the redshift distribution, appears poorly reproduced by ei- 
ther models. 

We then propose a new, unifled picture, in which the 
key feature is a link between the shape of the SEDs, in par- 
ticular the peak frequency of the synchrotron power distri- 
bution, and the bolometric luminosity. In particular, we pa- 
rameterize the shape of the SED in terms of its synchrotron 
peak frequency and assume a power law relation between 
peak frequency and bolometric luminosity. Adopting an es- 
timate of the bolometric luminosity function we can derive 
all the observable properties and again compare them with 
observations. 

We test all the scenarios against the best (whole) body 
of data now available: the radio and X-ray fluxes (uniformly 
measured by ROSAT, Urry etal. 1996) of the 1 Jy BL Lac 
sample (Stickel et al. 1991) and the radio and X-ray fluxes 
of the BL Lac sample derived from the Einstein Slew survey 
sample (Elvis et al. 1992; Perlman et al. 1996a). Based on the 
data and the discussion by Perlman et al., we treat the Slew 
survey sample available to date as a tentative, but quasi- 
complete one. 

In Section 2 we present the radio and X-ray samples of 

t We note that the definition of XBL and RBL by G&P is based 
on the X-ray flux in the 0.3—3.5 keV Einstein IPC band (in 
erg/cm^ /sec) and the radio flux at 5 GHz (in Jy). According 
to this difli'erent definition, the value of —5.5 adopted here corre- 
sponds to ~ —10.8 
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BL Lacs against which we tested the predictions of the mod- 
els. The two alternative "radio and X-ray leading" scenar- 
ios, as they will be called hereafter, are outlined in Section 
3, while we describe the simulation procedure adopted for 
the comparison with the observed samples in Section 4. In 
Section 5 we present and compare the predictions of the two 
models with observations and briefly discuss the effects of 
evolution. The assumptions and predictions of the new "uni- 
fied bolometric" model are described in Section 6. Section 7 
summarizes our conclusions. 

The values Ho = 50 km/sec/Mpc and go = 0.0 have 
been adopted throughout the paper. 



2 THE REFERENCE SAMPLES 
2.1 The 1 Jy sample 

The complete 1 Jy BL Lac sample was derived from the cat- 
alog of radio selected extragalactic sources with Fsghz > 1 
Jy (Kiihr et al. 1981) with additional requirements on ra^ 
dio flatness (ajj < 0.5, with F^, oc optical bright- 

ness (my < 20) and the absence of optical emission lines 
(EWa < 5 a, evaluated in the source rest frame) (Stickel 
et al. 1991). This yielded 34 sources matching the criteria, 
26 with a redshift determination and 4 with a lower limit 
on it (Stickel etal. 1994). It is the largest complete radio 
sample of BL Lacs compiled so far. 

We computed the luminosities using the monochromatic 
1 keV fluxes measured by ROSAT (Urry ct al. 1996, where 
the listed values are those from fits with Galactic absorp- 
tion) and the 5 GHz values from the Kiihr et al. catalogue. 
We considered only the subsamples of sources with at least a 
lower limit on the redshift. The fluxes were K-corrected us- 
ing a radio spectral index or = —0.27 and the average X-ray 
spectral index measured in the ROSAT hand, ax = 1.16. 



2.2 The Slew survey sample 

The Einstein Slew survey (Elvis et al. 1992) was derived from 
data taken with the IPC in between pointed observations. A 
catalog of 809 objects has been assembled with a detection 
threshold flxed at 5 photons. It does not reach high sensi- 
tivity, having a flux limit of ~ 5 x 10~^^ erg/cm^/sec in the 
IPC band (0.3 — 3.5 keV), but it covers a large fraction of 
the sky (~ 36600 deg^). Based on radio imaging and spec- 
troscopy Pcrlnian etal. (1996a) selected from this sample a 
set of 62 BL Lac objects (33 previously known and 29 new 
candidates) and, in a restricted region of the sky, a quasi- 
complete sample of 48 BL Lacs. The sources are almost com- 
pletely identified, and therefore constitute a quasi-complete 
sample of 48 BL Lacs. The redshift is known for 41 out of 48 
objects. This is the largest available X-ray selected sample 
of BL Lacs (Perlman et al. 1996a), with more than twice the 
number of sources contained in the fairly rich X-ray selected 
sample derived from the Einstein Medium Sensitivity Sur- 
vey (EMSS, Morris etal. 1991, Welter etal. 1994; Perlman 
et al. 1996b), extended to 23 objects by Wolter et al. (1994). 

K-corrections of the radio and X-ray monochromatic 
fluxes were computed using spectral indices ur = 0.0 and 
ax = 1.5, respectively. 



3 THE TWO SCENARIOS 
3.1 Different energy cut— off 

Giommi and Padovani introduced the idea that BL Lacs are 

a single population of objects whose SED can be character- 
ized phcnomenologically by the distribution of the values of 
the frequency at which the peak in the energy emitted per 
logarithmic bandwidth occurs (i.e. the peak in the vFv rep- 
resentation of the broad band energy distribution) for the 
putative synchrotron component. 

A good quantitative indicator of the overall broad band 
shape of SEDs is the ratio between X-ray and radio fluxes, 
conventionally taken at 1 keV and 5 GHz, equivalent to a 
characterization in terms of urx and oro (for a represen- 
tation of the relation between SED shape and for instance 
oro see Fig. 2 in Maraschi et al. 1995; Comastri et al. 1995 
tested also the orx — t'peak correlation). SEDs peaking at 
lower frequencies correspond to higher ratios (e.g. Maraschi 
et al. 1995). The X-ray to radio flux ratio is typically two 
orders of magnitude larger in BL Lac objects derived from 
X-ray surveys than in objects derived from radio surveys. 

It is clearly more convenient to work with flux ratios 
since they arc easier to determine then the peak frequency 
of a broad band energy distribution. One can even observa- 
tionally derive a distribution for the X-ray/radio luminosity 
ratio as if this is the relevant intrinsic quantity which char- 
acterizes the SED distribution. This is what we consider 
hereafter, following the G&P approach. 

Starting from these hypothesis, one can derive the sta- 
tistical properties of BL Lacs simply by assuming a radio 
(X-ray) luminosity function and a probability distribution, 
V, of the X-ray to radio luminosity ratio. In the follow- 
ing two subsections we present the assumptions on these 
two quantities, according to the radio leading model and its 
symmetric, in which the X-ray luminosity is the leading one. 



3.2 "Radio luminosity leading" 

The basic radio leading (G&P) hypothesis (Giommi & 
Padovani 1994; Padovani & Giommi 1995) is that radio 
selection would be objective with respect to the intrinsic 
spread of broad band spectral properties. The underlying 
idea is that the radio emission is only weakly affected by the 
properties of the synchrotron component at higher energies, 
such as the peak frequency. Therefore the radio selection is 
not expected to suflfer of any bias regarding the SED shape 
properties and sources with different X-ray/radio luminos- 
ity ratios are sampled from a common radio luminosity func- 
tion. 

In this approach, the X-ray counts (and X-ray lumi- 
nosity function) are easily predictable from the radio counts 
(and radio luminosity function) and the results can be com- 
pared with real data coming from surveys. It is worth re- 
minding that there are essentially no free parameters in the 
model. 

Since a radio selected sample would be unbiased with 
respect to the X-ray to radio luminosity ratio, G&P con- 
sidered the 1 Jy sample of BL Lacs (Stickel et al. 1991) and 
built the probability distribution V. To improve the statis- 
tics for the XBL-likc objects part of V, we considered, fol- 
lowing Padovani and Giommi (1995), the distribution de- 
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duced from the EMSS sample, weighted by a factor 1/10. In 
fact the ratio of XBL to RBL objects derived from various 
radio selected samples (1 Jy, Stickel etal. 1991; S4, Stickel 
& Kiihr 1994; S5, Stickel & Kiihr 1996) covers the range 
A^xbl/A^rbl = 1/16 — 1/6. Therefore, with the given hy- 
pothesis, this factor 1/10 represents the intrinsic ratio of 
XBL to RBL populations in the whole BL Lac class: in the 
radio leading scenario only a minority of sources would in- 
trinsically have the SED peaking in the UV range, while 
most of the energy distributions would peak in the ~ IR 
range. 

To calculate the luminosity function in the luminosity 
range io^*~^^ [erg/sec/Hz], G&P extrapolated the observed 
one to lower luminosities. This extrapolation is based on 
the beaming model scheme (see Urry etal. 1991), applied 
to the luminosity function derived from the 1 Jy sample 
data (covering the range 10^^ — 10^^ erg/sec/Hz, Stickel et 
al. 1991), in the non evolutionary case. We note that the 
predictions of the model are indeed strongly dependent on 
the number density of the low luminosity sources. 

Cleaxly, the assumption of a minimum radio luminos- 
ity LsGHz.min introduces a limiting luminosity in the X-ray 
band, Likev.min- However, it is worth noting that, because of 
the different spectral shape of XBL and RBL objects, their 
corresponding minimum luminosities will be different, with 

ilkeV,min|RBL < ilkeV,min |xBL ■ 

G&P found a good agreement between the predicted 
X-ray counts and those derived from all the main X-ray 
surveys (EMSS, EXOSAT High Galactic Latitude Survey, 
HEAO-1 surveys; see Fig. 3 of Padovani & Giommi 1995) 
and with the EMSS luminosity function. 

In their model, the observational evidence that XBL ob- 
jects are more numerous, apparently in contradiction with 
their assumption that XBL sources are intrinsically only 
1/10 of the entire population, is due to selection effects. In 
fact, given their -Fikev/-F5GHz ratios, XBL have a relatively 
low radio flux. Therefore, at the same X-ray flux level we 
are detecting together BL Lacs belonging to different parts 
of the radio luminosity function, with RBLs coming from a 
brighter but poorer part of it. 

Our aim here is to perform a further step in checking 
the radio leading scenario, i.e. to compute the distribution of 
radio and X-ray luminosities predicted for samples selected 
in the "leading" (radio) band also in the "secondary" (X- 
ray) spectral band, using the same inputs of G&P. 



3.3 "X— ray luminosity leading" 

The first alternative scenario that we tested is substantially 
"symmetric" to the radio leading one. It assumes that the 
common property of XBL and RBL is the X-ray luminos- 
ity. As already mentioned, this idea originally comes from 
the observational evidence that BL Lacs discovered in radio 
and X-ray surveys actually show similar X-ray luminosi- 
ties, while the radio luminosities typically differ by two 
three orders of magnitude (e.g. Maraschi et al. 1986; see 
also Fig. lb). This is in general true also if the distinction 
between the two types of sources is made in terms of the 
RBL/XBL classification as formally defined in Section 1 (see 
also Fig. la). 

Again, the starting points are the X-ray luminosity 



function and a probability distribution, V, of luminosity ra- 
tios, but in this case this should be observationally derived 
from an X-ray survey, that by construction is now supposed 
to be the objective one. 

An important assumption, which we had necessarily to 
take into account in the X-ray case, concerns the cosmologi- 
cal evolution of BL Lacs. While it has been so far established 
that RBLs show a slight positive evolution, consistent with 
no evolution at the 2cr level (e.g. Stickel et al. 1991; Wolter 
etal. 1994), recently Wolter etal. (1994) and Perlman etal. 
(1996b) confirmed the strong negative evolution in the X- 
ray band: X ray bright objects are much less luminous or 
common at high redshifts. In analogy with the radio case, we 
first evaluated the predictions of the X-ray leading model 
assuming no evolution. The results have been presented by 
Fossati etal. (1996) who showed that they are not compat- 
ible with the 1 Jy survey properties, predicting a number 
of RBL objects in large excess with respect to the observed 
one. 

Then, in order to be able to reproduce the observational 
results and in particular the redshift distributions (see be- 
low), here we assumed a negative luminosity evolution in 
the X-ray band. 

We derived the distribution V from the Slew survey 
X ray selected sample, but using only objects with redshift 
lower than 0.25, in order to better approximate a distribu- 
tion appropriate for a sample of sources at the same redshift. 

Unfortunately, the X-ray luminosity function for this 
sample is not (yet) available. For this reason we consid- 
ered the EMSS and its well studied X-ray luminosity func- 
tion (Morris etal. 1991; Wolter etal. 1994). We tentatively 
adopted an X-ray luminosity function matching that de- 
rived by Wolter et al. (1994), represented as a single power 
law with slope — —1.62. Since the normalization of the 
luminosity function would substantially affect only the to- 
tal number of sources, but neither their relative number nor 
the average luminosities, in order to recover the correct ab- 
solute number of objects we fixed the normalization (for 
the case with best results reported in Table 1) at the value 
log$o(iikcv) = —36.3. This value is within the confidence 
range given by Wolter etal. (1994), with Likev correspond- 
ing to -Lo.3-3.5kev = 10*° erg/scc. 

The luminosity evolution parameter is r = 0.142, ac- 
cording to the definition L{z) = L{0) exp{—tLB /t), where 
tLB is the look-back time in units of to (for go = 0.0, to = 
1/Ho). The luminosity range is 10^*"^° [erg/sec/Hz]. Also 
in this case we stress that RBLs and XBL reach different 
values of the minimum radio luminosity (see Section 3.2). 



4 THE SIMULATIONS 

Given these assumptions, we computed the predictions of 
both models. We used a Montecarlo technique to simulate 
the distribution of sources in space and luminosity. This 

method is very convenient because it allows us to store "sin- 
gle source" attributes and not only to compute sample inte- 
grated average properties. In other words, we can actually 
simulate a catalogue and from it compute the desired aver- 
age quantities. 

For the "radio leading" scenario we used a third degree 
polynomial fit to the luminosity function published in Fig. 1 
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of Padovani & Giommi (1995)[J. In all cases simulations have 
been performed within the redshift range — 2. 

In Fig. 2a, b the observed LikcV /^sghz distributions are 
shown, taking into account only objects with redshift esti- 
mate. The dashed histograms represent RBL objects, the 
dotted ones XBLs. Note that there is a superposition of the 
two classes due the fact that the actual distinction between 
RBL and XBL is based on the flux ratio (while here the 
luminosities are considered). In Fig. 2c the distribution V, 
used in the calculations for the radio leading scenario, is 
shown. 

The actual computation can be summarized in five 
steps: 

(i) "draw" a redshift and a luminosity from the luminos- 
ity function assumed as primary (radio for the radio 
leading, X-ray for the X-ray leading scenarios); 

(ii) the "secondary" luminosity (i.e. in the other spectral 
band) is deduced by choosing a luminosity ratio by 
means of its probability distribution V; 

(iii) radio and X-ray fluxes are computed for the simulated 
source taking into account the "inverse" K-correction^ 

(iv) check if the simulated source would be detected in a 
survey with a given flux limit and, if so, store its set of 
parameters. 

The detectability criteria take into account the sky cov- 
erage AQ that in the radio case is a step function of the 
monochromatic 5 GHz flux, while in the X-ray survey 
it is an increasing function of the flux integrated over 
the whole IPC band (0.3 — 3.5 keV). In the X-ray case 
we then simulated N surveys (typically 15 — 20) each 
with a different flux limit Fi and a corresponding Afli. 
The results of the N surveys are then summed up to 
construct the full X-ray survey. 

The flux dependent sky coverage for the Slew survey 
has been derived from the diagram "exposure time" vs. 
"percentage of sky" published in Elvis etal. (1992). A 
limiting flux corresponding to a given exposure time is 
deduced considering a detection limit of 5 IPC counts, 
and a conversion factor from IPC counts to flux of 
3.26 X 10"" erg/cmVsec (Elvis etal. 1992), valid for 
values of ax ~ 1-2 and Nh = 3 x 10^° cm~^. In fact, 
in order to convert to monochromatic flux we adopted 
an X-ray spectral index ax = 1-2, a flducial aver- 
age value for the class of sources that we are consid- 
ering. This gives the conversion relation: log -Fikev = 
log F[o.3_3.5keV] — 17.776. It is important to stress here 
that we use an approximated sky coverage by using 
the information currently available. An exact correction 
cannot be computed, but plausibly this would mostly af- 
fect the normalization of the assumed luminosity func- 
tions. 



t The polynomial fit is expressed as log <& = a + 
61ogL5GHz,28 +c(logL5GHz,28)^ +'^(logi5GHz,28)^ Gpc'^ L'^, 
where 1/50112,28 = 10~^*L5ghz erg/s/Hz, and a = —25.204616, 
b = -2.13762'71, c = -0.0455836 and d = -0.0059641 
§ The "inverse" refers to the fact that usually the K-correction 
is meant to convert observer-frame to source— frame quantities, 
while here we are using it in the other direction. 



(v) finally, a detected source is classified as XBL or RBL 
according to its -Fikcv/^sGHz ratio, as already defined. 



5 RESULTS AND COMPARISON OF THE 
MODELS 

The comparison of the two models with the 1 Jy and Slew 
survey samples are presented in Figs. 3 and 4, respectively. 
The average observed and predicted radio and X-ray lumi- 
nosities are reported in Table 1, together with the numbers 
of objects. The consistency of the predicted and observed 
distributions has been quantitatively estimated through the 
use of the Student's t-test for comparison of the means, 
and also of the more sophisticated survival analysis non- 
parametric and univariate methods (Feigelson & Nelson 
1985) allowing to take into account lower limits. Student's 
t-test has been performed on all average luminosities, while 
the survival analysis has been applied only to the cases with 
a significant number of objects, namely for the RBLs sub- 
sample in the "1 Jy" and XBLs in the Slew survey. The 
results of the statistical tests are reported as an entry in 
Table 1 for Student's t-test, and separately in Table 2 for 
the survival analysis. 

• Radio survey 

Let us consider first the comparison of the models with the 
1 Jy sample (Fig. 3). 

Both models predict luminosities for the XBL in both 
bands more than one order of magnitude larger than ob- 
served. However the extreme paucity of objects does not 
allow these discrepancies to be significant from a statistical 
point of view. The Student's t-test indicates that the radio 
leading scenario predicts a disagreement, at the > 94.0 per 
cent level, in the X-ray luminosity of the RBL, while the X- 
ray leading one has difficulties at reproducing the average 
luminosity of the XBL at the 91 per cent level. 

The relative number of objects detected in both scenar- 
ios is in reasonably good agreement with the observed one 
(5/29 and 3/27 instead of 2/32). 

•X-ray survey 

The predictions of the two models concerning the Slew sur- 
vey objects are instead quite different (Fig. 4). 

The radio leading scheme leads to incorrect X-ray and 
radio luminosities of XBL: from survival analysis 97.2 and 
99.96 per cent levels (Table 2), respectively, even worst ac- 
cording to the Student's t-test (Table 1). In particular, in 
the radio leading scenario the average radio luminosities of 
the RBL and XBL subsamples tends unavoidably to be the 
same: as immediate consequence of the "radio leading" hy- 
pothesis, all BL Lac objects should share the same range of 
radio luminosity, in contrast with observations of complete 
samples (see also Fig. lb). The average radio luminosity of 
XBL is also overestimated within the X-ray leading sce- 
nario (> 95.6 per cent level). Again, nothing can be said 
from these tests on the distributions of the less numerous 
source population. 

The predicted numbers of XBLs and RBLs in the Slew 
survey (43/7 and 43/7) are in good agreement with the ob- 
served one (40/8). 

Summarizing, in terms of source numbers and average lu- 
minosities the X-ray leading model is better than the radio 
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leading one, sharing with it the problem of XBL luminosities 
in radio surveys. Nevertheless, it should be noted that XBL 
objects in the 1 Jy survey are only 2 and so their measured 
properties could mis-represent the true ones. 

• Redshift distribution 

Finally, there is a significant disagreement on the red-shift 
distributions (see the panels of the first columns in Figs. 3 
and 4): both models predict a sharp peak of simulated RBL 
sources at low redshift, in contrast with the observed flat 
distributionsp|. 

This is a bmlt-m feature of this kind of approach in 
which we take the relative number of the the two types of 
sources fixed for any value of the leading luminosity. 

Note however that, as expected, the inclusion of the neg- 
ative evolution in the X-ray leading model tends to improve 
the agreement with observations. Indeed the redshift dis- 
tributions obtained in the X-ray leading non evolving case 
were qualitatively similar to the ones of the radio leading 
scenario (see Fossati et al. 1996). In fact in the 1 Jy sur- 
vey XBLs are sharply concentrated at lower redshifts, while 
RBLs show a "flatter" distribution. These trends are closer 
to the observed ones than those derived in the radio lead- 
ing scenario (Fig. 3). The same is true for the Slew survey 
(Fig. 4). 

The use of the survival analysis allowed us to perform a 
quantitative analysis. The comparison of the only RBLs in 
the 1 Jy survey indicates that the predicted distributions are 
not compatible with the observed one at the 95.80 and 93.45 
per cent level for the radio and X-ray leading scenarios, re- 
spectively. A much more significant disagreement is found 
in the Slew survey (considering only XBLs), which is incon- 
sistent with the predictions of the radio leading scenario at 
the 99.96 per cent level. 

5.1 Evolution 

This last point on the redshift distribution is particularly im- 
portant. In fact it indicates that evolutionary efi^ects should 
be taken into account. However it is possible that it is the 
distribution V which changes with redshift, rather than the 
luminosity in one of the two bands: clearly one can re- 
express the different evolution in the two bands, by saying 
that V evolves. 

From an observational point of view, if one considers dif- 
ferent intervals in redshifts in the Slew survey sample there 
are 5 RBL and 7 XBL objects at high redshifts, while these 
numbers become 3 and 26 at low redshifts (taking a dividing 
value at z = 0.25). The only two XBL sources with known 
redshift present in the 1 Jy sample are both at z < 0.1. 

We then considered an evolving probability distribution 
V, i.e. we introduced a redshift dependence on the relative 
ratio of the two populations of BL Lac objects. Basically V 
has a bimodal shape, described with two gaussians centered 
at log(I/ikGv/I'5GHz) = —6.5 and —3.5 for RBLs and XBLs, 

^ Note that the inconsistency between the distribution predicted 
by the radio leading scenario and the one observed in the 1 
Jy sample is in agreement with the assumption of no evolution 
adopted here. If we had introduced the sUght positive evolution 
apparent from "1 Jy" data we would have obtained a redshift 
distribution completely consistent with the observed one. 



respectively. The relative normalization of these gaussians 
evolves as NxBh/N-RBL oc (1 + z)~'^ (with 7 > 0), in such 
a way that the RBL part of the distribution becomes more 
and more dominant at higher redshifts. With this evolution 
of V we can induce the observed difference in the redshift 
distribution of XBLs and RBLs. 

In this case the luminosity function is not self- 
consistently derived. Instead, we have to deal with a tenta- 
tive luminosity function, in the sense that it is by no means 
consistent with the evolution that we are introducing. It is 
therefore necessary to derive "backwards" how the evolution 
of V translates in terms of luminosity function evolution. 

We considered only a limited region of the parameter 
space of 7 and A'xbl/A^rbl|{z=o)) without finding a solution 
which improves the results already obtained. Actually, in 
the case of the radio leading scenario, there is a hint that 
better solutions are obtained by changing the parameters 
progressively towards the no evolution case, as in the G&P 
model. 

While we cannot exclude that a suitable solution could 
be found following these assumptions, we decided to consider 
a new approach "a priori" , described in the following section. 



6 UNIFIED BOLOMETRIC APPROACH 
6.1 The idea 

Partly motivated by the difficulties found with the models 
discussed so far, we tried a completely different way of facing 
the problem. Furthermore, we wanted to take into account 
the properties of the overall observed SED of BL Lacs, as 
well as recent indication of a possible link, along a continuous 
sequence, between the SED shape and the source luminosity 
(e.g. Sambruna et al. 1996, Ghisellini et al., in preparation). 

The fundamental hypothesis of the model is still that 
XBL and RBL sources are manifestations of the same phys- 
ical phaenomenon. The new ingredient is the idea that XBLs 
and RBLs are different representatives of a spectral sequence 
that can be described in terms of a single parameter, which 
we identify with the bolometric luminosity of the synchrotron 
component, Lboi,sync. Both "flavours" of BL Lac objects 
share the same bolometric luminosity function and the SED 
properties depend strongly on this quantity. 

The main positive feature of this approach is that it 
offers a more direct interpretation in terms of the physical 
properties of these sources. In fact: 

a) the assumptions are largely independent of the details of 
the observed statistical samples and moreover they are 
not based on the choice of a leading spectral band; 

b) there is an immediate connection between the parameters 
of solutions "acceptable" from a statistical point of view 
and the physical conditions of the emitting plasma (see 
Section 7). 

The relation between the bolometric luminosity and the 
SED has been based on observed trends. More luminous ob- 
jects seem to have RBL-type spectral properties, with the 
peak of the energy distribution in the IR-optical range and 
Compton dominated soft X-ray spectra. 

The less luminous sources tend to display XBL-like 
SEDs, steep soft X-ray spectra and a synchrotron compo- 
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nent peaking in the UV-soft X-ray band (e.g. Fig. lb) . The 
distribution between this two extremes seems to be contin- 
uous and we propose that it is governed by the now leading 
parameter Lboi.sync. Therefore, if we characterize the SED 
with the frequency at which the (synchrotron) energy dis- 
tribution has a maximum, a fundamental (inverse) relation 
must exist between the bolometric luminosity and this fre- 
quency. The inverse dependence is qualitatively based on the 
trend found by Ghisellini et al. (in preparation). 

Interestingly, note that the different redshift distribu- 
tions of the two kinds of BL Lacs could be a natural outcome 
of this scenario. XBLs objects come from the lower (and 
richer) part of the luminosity function and so they would 
dominate at low redshift, but they would disappear at large 
distances despite the increase of the available volume. On 
the contrary, RBLs, even though coming from the poorer 
part of the luminosity function, would become predominant 
at higher redshifts, being still detectable. 

In the next Sections we discuss a parameterization of 
the observed SEDs, which allows us to quantify the predic- 
tions of this model. 



6.2 SEDs parameterization 

In order to reproduce the basic features of the observed 
SEDs we considered a simple two component model. The 
synchrotron radio to soft X-ray component is represented 
with a power law in the radio domain, with spectral index 
Ua, smoothly connecting, at z^junct, with a parabolic branch 
ranging up to fpeak- Beyond Vpenk the synchrotron compo- 
nent steepens parabolically. 

The hard X-ray Compton component is simply repre- 
sented with a single power law with spectral index a^. The 
normalization at 1 keV of this hard component (I/ikeV,Comp) 
is kept fixed relative to the synchrotron one at 5 GHz 
(isGHz). This immediately implies a correlation between the 
radio and X-ray luminosity in sources where the emission 
at 1 keV is dominated by the Compton emission. 

The parameters describing this representation of the 
SED are then five: Qs, ah, the hard (Compton) X-ray nor- 
malization K= log(Likev,Comp/I'5GHz), ^-junct and either the 
width of the parabolic branch a or Vpcak- This is because 
we require that the parabolic branch matches the flat radio 
component at t'junct also in its first derivative, and therefore 
we cannot fix independently both its width and peak fre- 
quency. A schematic representation of this parameterization 
is shown in Fig. 5, while its analytical expression is reported 
in Appendix A. 

For simplicity, we use the luminosity at the peak of the 
synchrotron component, I/boi,sync instead of the bolometric 
luminosity. In fact, for various complete samples of BL Lacs, 
Sambruna et al. (1996) find that the ratio between I/boi,sync 
and i-'peakipeak,sync is consistcnt with being constant, ~ 8. 

We then specify the relation between z^pcak (the fre- 
quency at which the peak in the i^Li, representation occurs) 
and the value of VpcakF' (i/pEak) {F' is the flux in the source 
frame). This is the key physical relation of the proposed 
model. For simplicity we considered a simple power law de- 



pendence l-poak = t'peak,o(ipoak/-£'poak,o) (»? > O)^- We 

then have three more parameters: Lpcak.o and Upcuk.o- 

In Fig. 6 a set of SEDs following the adopted pa- 
rameterization is shown. The dashed lines represent those 
with log(I/ikev/-'j5GHz) > —5.5, roughly corresponding to 
XBL-type objects, while the opposite is true for the solid 
ones, RBL-type objects. Remember that the distinction 
is actually defined in terms of the flux ratio, and that 

log(Fikcv/-F5GHz) = log(I/lkcv/i5GHz) + («« - Qx) log(l + 

z). 

To summarize: the adopted parameterization requires 8 
inputs, 5 describing the shape of the SED (qs, an, K, fjunct, 

t^peak 

) and three its relation with the luminosity (i^pcak,o, 
ipeak.o, ff)- However, the family of SEDs is completely de- 
termined by a set of 7 parameters, because the flrst 5 reduce 
to 4 + 1, (cfs, cch, K, t'junct, i^pcak(i^pcak,o, ipeak.o, ??))• lu ad- 
dition we must specify the normalization $o and the slope 
Q4> of the adopted luminosity function. 

6.3 SEDs: parameterization vs observations 

Despite the number of parameters, there is only a limited 
freedom in the choice of their values, which are constrained 
by observations. 

— Os only takes values in a narrow range around 0, and we 
perform our simulations with the fiducial values 0.0, 0.1, 
0.2. 

— The value of oh is constrained by the assumption that, 
in the case of extremely hard X-ray spectra, we are only 
looking at the Compton component. We then considered a 
single value of ah ~ 0.7 equal to the slope of the fiatter X- 
ray spectra of RBL (Lamer et al. 1996, Urry et al. 1996), 
for which we can think that we are measuring only their 
Compton emission. 

— On the same basis, the relative normalization between the 
radio and the X-ray components can be constrained by 
the value of log(I/ikcv/i5GHz) of extreme RBL type ob- 
jects, i.e. in the range [—8 ; —7]. We adopted the values 
K = -7.0 and K = -7.5. 

— 77 can be limited by the ranges of observed Lpcak and i/pcak, 
which imply values of rj between ~ 1 and 2. 

— i^junct is the parameter less determined by observations, be- 
ing difficult to identify a spectral break even in the case of 
the best sampled and simultaneous SED. We considered 
values between 10^" and 10^^ Hz. In the context of this 
parameterization the actual importance of fjunct is indi- 
rect: its position affects the width of the parabolic branch 
extending up to the X-ray band. For a fixed value of fpoak 
the higher is z/junct , the narrower will be the parabola (i.e. 
the steeper the cut-off). 

Let us consider now some observed properties of com- 
plete BL Lac samples, which allow us to check the goodness 
of our parametric representation. 

In Fig. 7, the classic aox vs ano plane for BL Lacs 
is shown, with points representative of the 1 Jy, Slew and 
EMSS samples. The line crossing the diagram is the locus of 



II We note that hereafter -Lpeak indicates the equivalent of a bolo- 
metric luminosity, i.e. uLi, evaluated at I'peak- 



8 G. Fossati, A. Celotti, G. Ghisellini, L. Maraschi 



points derived from our "best fit" SEDs. It basically repro- 
duces the observed pattern of two different branches, wliere 
RBLs and XBLs are concentrated, apparently with no over- 
lapping. The fact that the "synthetic path" does not extend 
beyond aox ^ 0.8 is related with the narrow range of pa- 
rameters (and tight relation among them) adopted in the 
samples simulation, namely to the minimum allowed bolo- 
metric luminosity. 

It is also possible to compare the SEDs with the ob- 
served average spectral energy distributions of RBLs and 
XBLs, taken from Sambruna et al. (1996). Their data points, 
given in flux units and converted in luminosities (consider- 
ing a typical redshift of 0.5 for the RBLs and 0.2 for XBLs), 
are plotted in Fig. 6 together with our SED family. Also in 
this case the agreement is satisfactory. 

Finally, also the LikeV vs I/sghz diagram (see Fig. lb) 
can be interpreted in terms of the SEDs of RBLs and XBLs. 
There is a fixed relation connecting the synchrotron radio 
luminosity and the inverse Compton X-ray luminosity and 
this can be better observed in RBLs, where more likely the 
Compton component dominates the X-ray emission. RBLs 
then occupy a locus following this correlation. The XBLs 
on the contrary lie away from this line because their X-ray 
emission is progressively dominated by a different compo- 
nent, the extreme part of the synchrotron spectrum. 

We indeed verified that the luminosity correlation for 
RBL seen in Fig. lb is not purely due to the common de- 
pendence on redshift of iikev and I/sghz- A Kendall's t 
test gives a correlation at >99.99 per cent level, which is 
still present, even though at the 96.8 per cent level, when 
the redshift dependence is excluded. Furthermore, we tried 
to establish the real nature of the lack of observed sources 
in the region around logL5GHz=33 and log LikoV=28.5. The 
results are shown in Fig. la, which is the equivalent plot, 
where fluxes (and not luminosities) are considered. The ver- 
tical and oblique lines represent the fiux limit of the 1 Jy 
survey and the ratio of the fluxes deflning XBL and RBL 
(i.e. —5.5), respectively. The grey area of the diagram shows 
the position in this flux-flux plane of putative sources with 
luminosities in the empty area of Fig. lb (indicated with 
similar shading) when we imagine tomove them in distance, 
covering a range of redshift between — 0.6. It can be seen 
that there is no apparent observational bias against detec- 
tion of sources with those luminosities: this indeed suggests 
that the lack of objects is plausibly linked to their intrinsic 
properties. This fact seems therefore to qualitatively agree 
with the shape of the SED assumed in the bolometric sce- 
nario. 

Clearly, none of the above checks concerns the num- 
ber density of objects. We adopted a luminosity function 
for I/pcak inspired to that calculated by Urry & Padovani 
(1995) (see their Fig. 13) for the 1 Jy and EMSS BL Lacs. 
Their (bivariate) luminosity function is obviously affected 
by selection effects both in the X-ray and radio bands, 
but constitutes the available distribution closest to a truly 
'bolometric' one. The normalization has been chosen in the 
range log ^(I/peak = 45.0) = -44.8 to —44.5 and the index 
a<s> = —1.90 to —2.3: also on these parameters we have a 
limited freedom. Note that the normalization is not actu- 
ally an interesting parameter and, as already mentioned in 
Section 3.3, varying $o only affects the absolute number of 
sources, neither the relative number of the two kinds nor the 



average luminosities. Finally, we stress that for simplicity no 
evolution has been assumed in the bolometric scenario. 

To conclude this Section, we can say that our schematic 
representation of the SEDs well reproduces the basic prop- 
erties of the observed BL Lacs broad band spectra. We have 
9 parameters and 14 observational quantities to match (8 
average luminosities and 6 source numbers). 

6.4 Results 

Even a simple "flrst attempt" set of the 9 input parameters 
gives a surprisingly good output. In Table 1 the best results, 
found after a systematic check for a grid of parameter values 
within the allowed ranges, are reported. The corresponding 
input parameters are listed in the note to the Table. 

The Student's t-test implies that the model predicts 
too high luminosities for the RBL both in the radio and X- 
ray band (at the 99.9 and 98.6 per cent level, respectively, 
see Table 1) when compared with the observed 1 Jy Survey 
distributions. The result is confirmed by the survival analy- 
sis (sec Table 2). For the Slew survey instead no significant 
disagreement has been found. 

Furthermore, a very positive consequence of this sce- 
nario is that it implies a qualitatively 'better' redshift dis- 
tribution, as can be seen in Figs. 3 and 4, with respect to the 
previous schemes. In particular we stress that the bolometric 
scenario is the only one which predicts, as direct result of the 
SED dependence on the bolometric luminosity, that the red- 
shift distributions of the two populations can be such that 
RBL would tend to dominate at high redsliifts, while the 
opposite would be true at low z (see Fig. 4). This behaviour 
could be seen as equivalent to either a redshift dependence 
of the X-ray-to-radio-fiux ratio distribution V or a negative 
X-ray luminosity evolution. 

There is still a problem with RBLs in the radio survey. 
Even though the shape of the z distribution is correctly flat 
(and this is the major difficulty of the two other models), it 
extends well beyond z = 1, which leads to a disagreement 
significant at the 92.7 per cent level (survival analysis). We 
believe this fact is directly connected with the excessive av- 
erage radio luminosities predicted by the model. 

In order to explore the consistency of the shape of the 
predicted and observed redshift distributions, we applied the 
survival analysis test only up to 2 = 1, and did not find in- 
consistency between the two. with the observed one. Clearly, 
this result only indicates that the shape of the predicted 
distribution can resemble the real one over the considered 
redshift interval. However, as we discuss in Section 7, the 
excess of objects predicted at higher redshift (and/or their 
overestimated radio power) could be understood in an even 
broader unifying picture. 



7 DISCUSSION AND CONCLUSIONS 

We have considered three different scenarios, which as- 
sume that BL Lac objects constitute a single population of 
sources, with different spectral energy distributions (mainly 
different peak frequencies of the synchrotron component), 
which originally mislead into thinking of a separation of BL 
Lacs into two classes. 
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A representation of the main features of these scenarios 
is shown in the cartoons of Fig. 8a,b,c. 

Two of these models (namely the radio and X-ray lead- 
ing ones; see Fig. 8a,b) are symmetric, and differ in the 
choice of the leading parameter, i.e. the radio and X-ray 
luminosities, respectively. The models are constructed from 
observed properties, like the luminosity functions and the 
probability distribution V of the flux ratios, with basically 
no free parameters. They can correctly predict a significant 
number of quantities in agreement with current observa- 
tions. However they (and in particular the radio leading one) 
fail to reproduce some of the observed distributions. 

We then proposed a unified bolometric model (see 
Fig. 8c), whose main characteristic is to link the bolometric 
luminosity with the energy of the synchrotron cut-off. This 
scenario is based too on observational trends, however it uses 
a schematic parameterization of the BL Lac SED as well as 
a semi-empirical luminosity function for the bolometric lu- 
minosities. We stress here that while it is true that the bolo- 
metric scenario contains a significant number of (quite con- 
strained) parameters with respect to the two other models, 
it also predicts the statistical and cosmological properties of 
the two 'types' of sources directly from a unified description 
of their spectral distribution. And, despite the rigid formu- 
lation of the one-to-one correspondence between the SED 
properties and the luminosity, in this scenario the main ob- 
servational data can be reproduced at least with the same 
'accuracy' as the other two models. Particularly interesting 
is the prediction of the different XBL and RBL redshift dis- 
tributions, which tend to favour the detection of objects of 
the first class at low redshifts and of RBL at higher z. 

The predictions of the bolometric model of luminosities 
of sources detected in radio surveys higher than observed is 
plausibly related to the problem with their redshift distribu- 
tion. We note, in fact, that the inclusion of Highly Polarized 
Quasars (HPQ) (i.e. emission line Blazars) in the LikcV- 
isGHz and colour-colour distributions shows an interesting 
continuity of their properties with those of BL Lacs (see 
Figs. 1 and 7, where stars represent HPQ). In particular, 
the significance of the correlation between Lsghz and iikcv 
(see Section 6.3) increases to 99.99 per cent (already exclud- 
ing the common dependence on redshift) when also the HPQ 
of the Impey & Tapia (1990) sample are included. The mor- 
phology of the extended radio emission of some RBL (see 
e.g. KoUgaard etal. 1992; Perlman & Stocke 1993), and the 
radio luminosity functions of RBLs and FSRQ (a class of 
Blazars substantially equivalent to HPQ) also show continu- 
ity (Maraschi & Rovetti 1994). All these pieces of evidence 
lead to the suggestion that there is actually a remarkable 
progression in properties, in a luminosity sequence XBL- 
RBL-HPQ, which is also a sequence of increasing impor- 
tance of emission lines. If a fraction of the sources with the 
highest bolometric luminosities could be indeed "classified" 
as HPQ rather than RBL (in a way analogous to the transi- 
tion from XBLs to RBLs), this fact could plausibly explain 
the discrepancies in the high luminosities and redshift dis- 
tributions predicted by the model. 

As already stressed, a very interesting aspect of the 
"bolometric" approach is that it is based on the relative 
dependence of two quantities which are strongly related to 
the physical properties of the emitting plasma, namely the 
emitted luminosity and the cut-off in the synchrotron spec- 



trum. One can therefore speculate on the physical origin of 
this dependence. 

Probably the simplest interpretation of this link is that 
it is directly related to the particle cooling. The more radia- 
tion is emitted the more particles loose energy, with con- 
sequent decrease in the cut-off frequency (which refiects 
the maximum energy of the emitting particles). If the en- 
ergy of the particles emitting at the peak of the fiux distri- 
bution is determined by the escaping time from the emit- 
ting region, this can be estimated by equating the cooling 
and a typical escape timescales. This leads to the relation 
J^poak oc (Lboi/i?)^^, where R is the size of the emitting re- 
gion. This dependence is reasonably close to the one which 
predicts the best results of the model (r; — 1.5). Clearly this 
interpretation is not unique and one can as easily envisage 
physical scenarios where a different relation between these 
two quantities is expected (e.g. in the case the source size is 
proportional to the central object mass). 

Furthermore the i'peak(iboi,syn) relation can be also at- 
tributed, in the context of inhomogeneous jet models, to a 
possible dependence of the length of the coUimated jet on 
the radiated power: more powerful jets could extend farther 
and originate outer radio regions responsible for the X-ray 
Compton component. 

To conclude, we consider here the predictions of the 
models examined on the results of deeper surveys. In par- 
ticular we compute the expected ratio of objects of the two 
"flavours" in more sensitive radio and X-ray surveys, as a 
function of their fiux limits. The results are reported in Ta- 
ble 3 |**| , as the relative ratio A^'rbl /Nxbl normalized to the 
value predicted for the 1 Jy and Slew surveys. 

Note that this ratio is a function of the flux limit even 
in the radio survey for the radio leading scenario (and X-ray 
survey for the X-ray leading one) . This is due to the different 
dominant range of redshifts at different flux limits, which 
causes objects at the border of the definition of XBL/RBL 
to "move" into the RBL class, because of the term {an — 
ax) log(l -I- z) (always < 0) in the relation between flux and 
luminosity ratios. 

Another factor to take into account in the interpreta- 
tion of these predictions is the influence of the minimum 
luminosities. As mentioned in Section 3.2, a minimum lu- 
minosity in the "leading" band leads to different minimum 
luminosities for RBL and XBL in the "secondary" band, due 
to their intrinsically different SEDs. Therefore, at flux limits 
which allow the detection of the faintest sources, only ob- 
jects with SEDs which favour very low fluxes in that band 
will be detected. This effect is observable mainly in radio 
surveys in the X-ray leading scenario and in X-ray surveys 
in the radio leading one. 

From Table 3, one can clearly see the opposite trends 
expected from the radio and X-ray leading scenarios in radio 
surveys. Distinct features of the bolometric scenario are a 



** Browne & Marcha (1993) pointed out the possible effect of 
misidentification of faint BL Lacs caused by the lack of contrast 
of the active nucleus with the host galaxy. It should be noted 
here that, while we were able to ignore the Browne & Marcha 
selection effect for the results of the Slew survey, this would be 
progressively relevant in surveys with decreasing flux limit, as the 
ones we are simulating here. 
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rapid decline in the fraction of RBL in radio surveys and a 

weak decrease in X-ray surveys. This behaviour follows from 
the one-way relation between RBL and XBL SED properties 
with luminosity. 
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8 APPENDIX A: SED PARAMETERIZATION 

The analytical expression of the parameterized SEDs in the 
"bolometric" model is given by: 



= I3{x - 9.698) + ipibGHz) 

i>s,2{x) = log(i^Fi,)s,2 = 

= - ((a; - a;peak)/f)^ + 0.25it^/3^ -|- tps,i{x) 

iph{x) = \og{vF^)h = 

= (1-Qh)(a;- 17.383) +Vcomp(lA:ey) = 

= (1 - ah){x - 17.383) + ipiSGHz) + K + 7.685 



where x = logz/, and /3 = (1 — as), and cr = 
(2(a::pcak - a;junct)//3)^''^. tps,iix) and -<l)s,2{x) are defined on 
the frequency ranges x < xjunct and x > xjunct, respectively. 

The SED expression can be re-scaled as a function of 
the value of Lpeak = lo^'^^p"'-''^, and in particular the lumi- 
nosities at 5 GHz and 1 keV can then be written as: 



tpi5GHz) = 

V'sync(lfcey) = 

V'Comp(lA:e\/) = 

V'tot(lfceV) = 



(1 -I- 0.5?7/3)V'peak - 0.5/3(Xpeak,0 + 
+?7V'peak,0 - Xjunct) - PiXjunct - 9.698) 
^(aJpeak) - ((17.383 - Xpeak)/lT)^ 
i>synri5GHz) + K + 7.685 

log(10'^=^" -I- 10*°°°p) 
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TABLE 1 
Summary of R,f,stit,ts^ 



1 Jy Survey 



SL 



ew survey 



RBL 



XBL 



RBL 



XBL 





Ntot 


Nrbl 


NxBL 


LbGHz 


LlkeV 


LbGHz 


LlkeV 


Ntot 


Nrbl 


NxBL 


LbGHz 


LlkeV 


LbGHz 


LlkeV 


Observed 


34.0 


32.0 


2.0 


34.05 
(0.76) 


27.30 
(0.30) 


32.10 
(0.30) 


26.84 
(0.23) 


48.0 


8.0 


40.0 


33.74 
(0.88) 


27.56 
(0.96) 


31.60 
(0.49) 


27.53 
(0.64) 


Radio leading 


29.1 


26.7 


2.4 


33.73 
(0.97) 
[ ••• ] 


26.83 
(1.24) 
[ 94.50 %] 


33.54 
(0.86) 
[ ••• ] 


28.93 
(1.05) 
[ ••• ] 


50.7 


7.1 


43.6 


33.10 
(0.99) 
[ ••• ] 


27.31 
(1.06) 
[ ••• ] 


32.43 
(0.98) 
[ 99.99 %] 


28.20 
(1.01) 
[ 99.83 %] 


X— ray leading 


33.7 


29.0 


4.7 


34.13 
(0.80) 
[ ••• ] 


27.51 
(0.83) 
[ ••• ] 


33.34 
(0.78) 
[ ••• ] 


28.50 
(0.78) 
[ 91.00 %] 


49.4 


6.7 


42.7 


34.10 
(1.03) 
[ ••• ] 


27.74 
(0.90) 
[ ••• ] 


32.07 
(1.10) 
[ 97.74 %] 


27.72 
(0.96) 
[ ••• ] 


Bolometric^ 


35.1 


33.3 


1.8 


34.87 
(0.78) 


27.83 
(0.33) 


32.43 
(0.46) 


27.31 
(0.21) 


47.2 


8.2 


39.0 


.33.57 
(0.76) 


27.64 
(0.21) 


31.85 
(0.49) 


27.28 
(0.35) 



; 99.93 %] [ 98.58 %] 



Poisson error should be considered on the source numbers. The variances on the average luminosities and the significance of rejection of the "null'* hypothesis as evaluated with the Student's 
t-test are reported in "( )" and "[ ]'', respectively. Dots are reported when the significance is lower than 90 

^The input parameters corresponding to this results are: = 0.2, aj^ = 0.7, K = —7.5, i^junct = 5 X 10-^'-' Hz, log i^peak,o = 16.5, log Lpeak,o = 44.5, r] = 1.5, <^o = —44.76, = —2.25 
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CM 



TABLE 2 
Survival Analysis Results^ 





Radio 
leading 


X-ray 
leading 


Bolometric 


1 Jy Survey^ 


z 


95.80 


93.45 


92.70 


LbGHz 






99.85 


LlkeV 


94.00 




94.80 


Slew Survey^ 


z 


99.96 






LbGHz 


99.96 


95.60 




LlkeV 


97.20 







^ A variety of statistical test has been applied: Gehan Gen- 
eralized with both permutation and hypergeometric variance, 
Peto &L Peto and Peto &l Prentice, all of them as implemented 
in the ASURV package (Rev. 1.1, La Valley, Isobe and Feigel- 
son 1992). The reported values are the simple averages of the 
results of these four tests, expressed in percentual probabil- 
ity. Dots are reported when the significance is lower than 90 

%. 

^In the "1 Jy survey" case the tests have been performed 
only on the RBLs subsample. 

'^In the "Slew survey" case the tests have been performed 
only on the XBLs subsample. 



TABLE 3 
Model Predictions^ 



Radio Survey X— ray Survey 



Flux limit^ Radio X— ray Bolometric Flux limif^ Radio X— ray Bolometric 

(FsGHz) leading leading leading leading 



1.00 1.0 1.0 1.0 1/1 1.0 1.0 1.0 

0.10 1.3 0.89 0.29 l/lO 1.63 1.13 0.95 

0.01 1.5 0.62 0.09 1/100 2.43 1.26 0.88 



^The numbers represent the predicted Nr^bl /Nxbl ratio normalized to the value predicted for the 1 Jy or Slew 
survey 

^ Radio fluxes are expressed in Janskys 

'^The flux limit is expressed relative to the Slew survey one. The sky coverage used is simply that of the Slew 
survey shifted, in logarithm, by the appropriate factor 
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Figure 1. (a) The log i^ikcV vs log Fsghz plane for different samples: 1 Jy, EMSS, Slew survey, and the complete HPQ sample from Impey 
&; Tapia (1990) (all the corresponding symbols are defined in the inlet in panel (b)). The vertical line represents the IJy flux detection 
limit, while the dashed line defines the separation between XBL and RBL. Finally, grey area shows the loci of points that potential 
sources shown in Fig. lb in the similar grey region, would occupy in a i^ikeV ^5GHz diagram for different redshifts (0 < z < 0.6). This 
shows that no obvious selection effects would avoid the detection of these sources, (b) Distribution of luminosities for the same source 
samples. The oblique line indicates the values of iikcv/^SGHz roughly corresponding to the definition of RBL and XBL. 
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Figure 2. (a,b) Distribution of the observed ratio LikevZ-^'BGHz for the two surveys. Dashed areas represent RBL, while dotted areas 
indicate XBL sources. The distribution ■P(-J/ikev/^BGHz) used in the computations of the radio leading scenario is shown in (c). It is 
obtained summing up that observed for the RBLs in the IJy sample and that of the EMSS sample, weighted by a factor 1/10. 
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Figure 3. Results of the simulations of the 1 Jy survey compared with observations. Prom top to bottom: observations, "radio leading" , 

"X-ray leading" and "bolometric" models, respectively. From left to right: histograms of the redshift distributions, radio and X-ray 
luminosities for the two classes. Again dashed and dotted areas represent RBL and XBL, respectively. The error bars are derived from 
the simulations. 
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Figure 4. Same as Fig. 3 for the simulated Slew survey. 
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^IkeV 
^5GHz 



5 GHz Vjunct Vpeak 1 keV 

Figure 5. Schematic representation of the parameterization of the SEDs, according to the "unified bolometric" model. All the basic 
parameters are shown. Note that the difference between radio and X-ray luminosities in the fLi^ representation differs from the relative 
normalization of the two component, defined as K, by the term log(i^ikev/'^6GHz)- 
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8 9 10 11 12 13 14 15 16 1? 18 19 20 

log u 

Figure 6. Set of SEDs used in the simulations of the "unified boloinctric" scenario. The dashed SEDs are those with log(Li]jeV / ^5GHz) > 
—5.5, corresponding to XBL objects. The sohd SEDs, on the contrary, correspond to RBL objects. Note that this subdivision is not 
strictly equivalent to the classification defined In terms of the fiux ratio, that differs by a term [{an — ax) log(l + z)] ~ —0.3. The squares 
Indicate the average observed SEDs of XBL and RBL, as reported by Sambruna et al. (1996). 
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Figure 7. Distribution of sources in the colour— colour diagram, i.e. in the apx vs a^o plane. Different completesamples are considered: 
EMSS, Slew and 1 Jy BL Lacs, and High Polarization Quasars (HPQ) from the Brunner et al. (1996) sample. All the corresponding 
symbols are defined in the inlet. The crosses connected by a dotted line show the spectral indices predicted by the adopted "best fit" 
parameterization of the SEDs. 
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(a) 




(b) 





Figure 8. Cartoons resuming the basic features and differences between (a) the radio leading, (b) the X-ray leading and (c) the new 
unified bolometric scenarios. The plane defined by the two axis shows the leading luminosity function. As shown in the panels in the 
first two schemes the SED distributions and populations ratios are independent from luminosity, while the bolometric approach links 
the type of source (i.e. the SED) with the luminosity. The schematic SED represented with thick and thin lines always refer to RBL and 
XBL, respectively, and the number of SEDs in each panel is suggestive of the relative number density of sources. 



